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Abstract 

The violent relaxation and the metastable states of the Hamiltonian Mean-Field 
model, a paradigmatic system of long-range interactions, is studied using a 
Hamiltonian formalism. Rigorous results are derived algebraically for the time 
evolution of selected macroscopic observables, e.g., the global magnetization. 
The high and low energy limits are investigated and the analytical predictions 
are compared with direct A^-body simulations. The method we use enables us to 
re-interpret the out-of-equilibrium phase transition separating magnetized and 
(almost) unmagnetized regimes. 

Key words: Vlasov equation, Hamiltonian systems, out-of-equilibrium phase 
transition 
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1. Introduction 

Systems with long-range interactions exhibit a fascinating feature of 

metastability: Starting from out-of-equilibrium initial conditions, the system 
violently relaxes toward a metastable state, often called Quasi-Stationary State 
(QSS). In this regime, macroscopic quantities reach values which substantially 
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differ from the corresponding thermodynamic equihbrium configuration. Al- 
though the QSS are only transient regimes, their lifetime have been shown to 
diverge with the number of bodies in interaction 3] . For this reason they pos- 
sibly correspond to the solely accessible experimental regimes. 

We consider a paradigmatic system with long-range interactions, the Hamil- 
tonian Mean-Field (HMF) ^ where particles on a circle are collectively inter- 
acting through a cosine-like mean-field potential. After a fast relaxation, the 
system typically enters a metastable regime in which the particles either aggre- 
gate into a large cluster (magnetized phase), or they spread almost homoge- 
neously around the circle (unmagnetized or homogeneous phase). In particular, 
an out-of-equilibrium phase transition between these two states occurs when the 
parameters of the initial conditions are varied Q . 

In this paper, we focus on both the violent relaxation process and the sub- 
sequent QSS regime. We use an algebraic framework based on a Hamiltonian 
formulation of the Vlasov equation for the HMF model. This Vlasov equation 
rules the evolution of the single particle distribution function in phase space 
(as a kinetic equation) and naturally arises when investigating the continuous 
version of the HMF model. As in the limit of infinite number of particles the 
system gets permanently frozen in the QSS phase, it is customarily believed 
that QSS can be interpreted as equilibria of the Vlasov equation. We exploit a 
Hamiltonian formalism of this Vlasov equation to derive analytical expressions 
for the global magnetization as function of time. This magnetization measures 
the aggregation of the particles on the circle. It is a macroscopic observable 
which is directly influenced by the microscopic, single particle trajectory. It is 
in general particularly cumbersome to bridge the gap between the microscopic 
realm of the many-body interacting constituents and the macroscopic world of 
collective dynamics. 

Using an expansion provided by the Hamiltonian framework, we here obtain 
rigorous results on the time expansion of relevant observables. These results are 
compared with direct numerical simulation. We consider in particular the high 
and low energy regimes which allow some simplifications in the expansions. In 
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addition, we characterize the aforementioned out-of-equilibrium phase transition 
which occurs in an intermediate energy range. This is achieved by monitoring 
the initial relaxation of the magnetization, as a function of relevant parame- 
ters of the initial distribution. The parameter space is hence partitioned into 
two regions, depending on the magnetization amount, a result which positively 
correlates with direct numerics 1^. 

The paper is organized as follows: In Sec. [5] we will review the discrete HMF 
model, presents its continuous counterpart and discuss the basic of the bracket 
expansion method. Section [3] is devoted to the presentation of the analytical 
results, with special emphasis to the high and low energy regimes. The out- 
of-equilibrium phase transition issue is also addressed. Comparison with direct 
simulations is provided to substantiate the accuracy of our predictions. 

2. Model and methods 

2.1. Lie-Poisson structure of the Vlasov equation 

We consider N particles interacting on a circle with the following Hamilto- 
nian: 



where {9i,pi) are canonically conjugate variables which means that the Poisson 
bracket giving the dynamics (Hamilton's equations) is given by 



In the continuous limit, we consider an Eulerian description of the system which 
gives the dynamical evolution of the distribution of particles f{6,p;t) in phase 
space via the following Vlasov equation: 



where V[f]{0) = 1 - M^[f]cose - My[f]sme. The magnetization M[f] = 
Mx + iMy is defined as 




(1) 





Of dV[f]df 
80 dO dp' 



(2) 




(3) 
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where the integrals are taken over [— tt, tt] x R. Equation ^ can be cast into a 
Hamihonian form where the (infinite dimensional) phase space is composed of 
the functions f{0,p) of ] — tt, tt] x K.. The Hamiltonian is given by 

and the associated Lie-Poisson bracket by 

rx ff ff d 6F d 5G d SF d SG\ 

^^'""^ = J J \d-pJfdeJf - deJId-pJin^^ 

for F and G two observables (that is, functionals of /). The functional deriva- 
tives SF/Sf are computed following the expansion : 

F[f + VP] - F[f] ^-^^dOdp + 0{^''). 

The Poisson bracket (O satisfies several properties: bilinearity, Leibniz rule and 
Jacobi identity (for more details, see Refs. 0, 0]). Its Casimir invariants are 
given by 

C[f] =11 c{f)d0dp, 



where c(/) is any function of f{9,p). In particular, the total distribution 
JJ fdOdp is one of such Casimir invariants and hence is conserved by the flow. 
The evolution of any observable F[f] is then given by 

F = {H,F}. (6) 

For instance, for F[f] = f{9,p), we recover Eq. Another convenient ob- 
servable to study is the magnetization M[f] given by Eq. ([3]): It quantifies the 
spatial aggregation of the particles. At low energies, the magnetization typically 
relaxes until it reaches an out-of-equilibrium plateau, around which it fluctuates 
(see Fig. [T]). In this case, the particles are trapped into the large resonance cre- 
ated by the finite magnetization, hence the name "magnetized state" (see upper 
panel of Fig. [2]). At high energies, the magnetization falls and fluctuates around 
zero (see Fig. [Ij , which means that the particles failed to organize collectively. 
This is called the "homogeneous phase" (see lower panel of Fig.[2|). 



4 



0.8 

0.6 

§ 0.4 
OS 

0.2 


-0.2 







Low-energy regime 




High-energy regime 

10 15 20 

t 



Figure 1 : Real part of the magnetization given by Eq. (|3]l as a function of time obtained by 
integrating the dynamics given by Eq. ^ for Mq = 0.6. The system reaches either a finite- 
magnetization for low energies {U = 0.4, in blue), or a low-magnetization for high energies 
{U = 3, in red). The plain lines refer to Af-body simulations (with A'^ = 10000), while the 
dotted lines come from the predictions given by Eq. lO for ko = 20. 

The dynamics given by Eq. ([6]) is deduced from the hnear operator Ti.. From 
the evaluation of the functional derivative of H with respect to / 

— = ^ - [/] cos e - My [/] sin 

we get the expression of TL : 

■H = {H,.} 



dOdpf p 



d Me- 



-ie 



, -■ (7) 

86 2i dpj Sf ^ ^ 

In the algebraic computations that follows, we make an explicit use of the 
linearity of Ti and Leibniz rule: 

n{F + aG) ^nF + aUG, 
n{FG) = FUG + {nF)G. 

N-body simulations (Lagrangian point of view): In order to compare the 
algebraic results with numerical ones, we integrate Eq. ([6|) via A'^-body simu- 
lations, which are obtained by considering a Klimontovitch Q distribution of 
particles 



fiO,p;t) 



1 ^ 
N ^ 



{e-e,{t))5{p-pm. 
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Figure 2: Phase space portrait of the system ((TJ once saturation has been reached for Mq = 
0.6: The low energy regime {U = 0.5, upper panel) is characterized by one large cluster of 
particles, whereas for higher energies [U = 1.7, lower panel), phase space is quite homogeneous, 
except for two clusters at it 2.2 (i.e. moving in opposite directions). 
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whose dynamics is ultimately reduced to Hamiltonian ([T]). Such simulations are 
used with a large number of particles (typically N — 10^) such that the interme- 
diate regime experienced by the A^-body simulations is close to the behavior of 
the Vlasov equations (at least for some observables like the magnetization) [sl- 

2.2. Bracket method 

The evolution of a selected observable F[f] given by Eq. ([6]) is obtained 
formally from the operator Ti. as 

oo j-k 

F[/](t) = e*«f^[/o]^^^H'=^[/o], (8) 

fe=0 

where fo is the initial distribution. Here we compute a finite number of terms in 
this series in order to obtain a Taylor expansion for the solution of the dynamics 
of F: 

kg J, 

fe=0 

where kg is the truncation parameter. Of course, this approximation is accurate 
up to some time t depending on /cq- A convergence over longer times is expected 
for increasing /co- Furthermore, the coefficients of the series H'^F are obtained 
recursively by applying Ti on the previous term Ti'^^^F. Finally, note that 
Eq. ([9]) yields an explicit dependence on the initial conditions, and it is not 
restricted to close-to-equilibrium initial conditions, thus being a useful tool to 
investigate the far-from-equilibrium violent relaxation of the system. 

We consider the subspace of functions composed by sums and products of 
the following elements (which are also functionals of /): 

buAf] = JJ d9dpfe"''p"', 

where (n, m) e Z x N. We notice that the main observables of the system such 
as the n-th order magnetization M„ = J J e^^^fdOdp or the momenta Pm = 
// P^^fdOdp of the system belong to this family. Furthermore this family is 
stable by the action of H given by 

Tn 

'Hbn,m — inbn,m+l + -^TT (^l.O^ri- l,m- 1 — l,O^Ji+l,m- 1 ) • (10) 



7 



We notice that only positive values of m arc involved in the iterations of the 
recursion relation since the second term is proportional to m. Taking into 
account the linearity and the Leibniz rule for Ti. mentioned in the previous 
section, the derivation of the short-time evolution ([5]) of a given observable F 
is computed algebraically as a sum of products of elements &„.,„. For instance, 
the magnetization is given by M — 61.0 and its first order evolution is obtained 
from Eq. PH)) : 



Of course, a satisfying approximation of the time evolution of any observable 
needs a large number of terms in the expansion ([9|). At a given time t, the 
number of terms necessary to obtain a reasonably good approximation of the 
dynamics depends on the initial distribution /o as it is shown in Fig. [T] where, 
at low energies, the accuracy extends to longer times than at high energies. 
In addition, we need to specify the initial distribution which will be used to 
compute bn,m{0) necessary to complete the computation of the approximate 
evolution. In the following sections, we use a waterbag distribution as initial 
condition. 

2.3. Initial conditions 

The waterbag initial distribution is a uniform distribution over a rectangle 
in phase space corresponding to the points {6,p) G [— A6', A^] x [— Ap, Ap]. The 
distribution fo{0,p) is equal to l/(4A6IAp) if {9,p) G [-AO, AO] x [-Ap,Ap] 
and zero otherwise. The values of bn,m at i = can be computed explicitly in 
this case and are equal to 



M{t) 




bn.m{0) 



{Ap} 



m+l -(-Ap) 



m+1 



sinc(nA0), 



2(m + l)Ap 
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where sinc(-) = sin(-)/(-). In particular, we notice that &n,m(0) = for m odd. 
The waterbag is characterized by two parameters (A0, Ap). Instead we consider 
the initial magnetization Mq and the energy U to label the initial conditions: 

Mi,si>i,o(0) = smc(Ae), 

In the following, we investigate the high energy t/ 3> 1 and low energy U 1 
limits for the initial distribution. 

3. Analytical results 

The first terms of the expansion for the magnetization M{t) given by Eq. ([9]) 
(for a waterbag initial distribution) are listed in Tab.[T]up to sixth order in time. 
We notice that the number of terms in the expansion increases exponentially, 
making such expressions difficult to handle in practice. In Fig.[Tl we notice that 
even with ko — 20 which involves approximately one thousand terms, a good 
agreement is observed only up to f = 2. In Fig.[3l the algebraic expressions for 
the magnetization obtained by Eq. ^ are plotted at different orders. Other 
than the initial regime, if one is interested in the intermediate regimes, the only 
hope is to find the governing rules behind this algebraic computations in order 
to draw some conclusions. This is the case for the low and high energy limits 
where the leading terms of the expansion can be extracted to all orders. These 
simplifications allow us to derive some dynamical properties of the system. 

Table 1: First terms in the expansion of the magnetization M{t) 
given by Eq. ^ for the waterbag initial distribution. 



iV2! 




(1 - sinc(2A6'))A/o/2 




Ap2 


-A/o/3 


i74! 


ApO 


(1-2 sinc(2A6') + sinc(2A6')2 - iM^ + AMq sinc(3A6l))A/o/4 




Ap2 


-2(3 sinc(2A0) + l)Afo/3 




Ap4 


Mo/5 
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Figure 3: Magnetization M{t) versus time for Mq = 0.2 and U = 0.6 obtained from Af-body 
simulations (dotted black curve) and using the algebraic expansions at various orders from 
the 4th to the 20th order. 







(1 + 3 sinc(2A6') + 64:M^ + 3 sinc(2A6l)2 + 26Mo sinc(3A6') 
+98M2 sinc(2A6l) - sinc(2A6')3 - SiM^ sinc(4A6l) 
-26Mo sinc(2A6') sinc(3A6l))Mo/8 




Ap2 


(-202 sinc(2A6')Mo - 51 sinc(2A6')2 + 58 sinc(2Ae') + 138Mo2 
-7)Mo/12 




Ap4 


(-239 sinc(2A6') + 23) Afo/30 




Ap*^ 


-Mo/7 



3.1. High-energy limit 

First we consider the high energy hmit, which corresponds to Ap ^ 1 in 
the initial waterbag. In this regime, since the kinetic term is dominant, the 
dynamics is driven by the reduced Liouville operator, which takes into account 
only the kinetic term 

nHE = II dddpfp^jj. (11) 
From Eq. pTI)l . the successive actions of H on bn,m is given by 
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Figure 4: Real part of the magnetization M{t) versus time for Mq = 0.6 and U = 20. The 
dashed red hne refers to direct simulations, while the solid blue one stands for the approximate 
solution l |12ll . 

For the waterbag initial distribution, it is straightforward to deduce the evolu- 
tion of the magnetization of order n: 

Mn{t) = M.n{0)smc{nApt). (12) 

The magnetization envelop exhibits a slow decay (as l/{Apt)) towards the 
asymptotic (equilibrium) state M — (see Fig. SJ. 

The profile obtained from A'^-body simulations is correctly interpolated over 
a finite time window by Eq. (|12p . As U is increased, the agreement gets better, 
even if deviations from Eq. (jl2p are observed at later times. Such a discrepancy 
is due to the cumulative effects of the neglected contributions in Ap (see Tab.[T]). 

It was reported in Ref. Q that for large values of the energy, and for any given 
initial magnetization, two large resonances spontaneously develop and effectively 
divide the available phase space into independent regions. Such resonances 
move in opposite directions, over the unit circle. Their velocity pr is identical 
in modulus and tends to grow as the energy is increased. The magnetization 
M{t) is mostly influenced by the instantaneous positions of the resonances. A 
snapshot of the positions of the particles obtained using A^-body simulations 
is depicted in Fig. [2] (lower panel). It reveals the two resonances moving in 
opposite directions (with velocity pr). The maxima of M{t) are obtained when 
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Figure 5: The velocity of the resonances (in the high energy regime) as a function of the 
energy U for Mq = 0.7. The circles refer to the velocities obtained numerically using A'^-body 
simulations. The vertical bars delimit the width of the resonances in the p-direction. The 
solid line is given by Eq. JTS}, and the dashed line by Eq. H14| l. 

the two resonances are aligned since the bunching of particles is maximum in this 
case. During two successive maxima of M{t) each cluster travels on a segment 
of length 2tt in 9, which takes a time 2n/pr- On the other hand, according to 
Eq. (|12p. two successive bumps in the magnetization are separated by a time 
interval 27r/Ap. This leads to 



which applies in the high-energy limit. The above prediction is compared with 
A^-body simulations in Fig. ([5]): The velocity (and corresponding width) of the 
resonances is plotted for different energies U (circles). The solid line refers to 
the analytical expression (|13p. We notice a very good agreement between the 
numerics and the prediction p3)) . As expected, as U decreases, some discrepancy 
is observed since the system approaches the phase transition. 

The above conclusion and in particular Eq. (jl2p can be also recovered using 
the following argument: In the high energy limit, the particles move essentially 
freely. The potential energy accommodates for just a small fraction of the total 
energy. Under this hypothesis, the individual phase 6 evolves as: 




(13) 



e{t) = 00 + Pot 
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where the index refers to the initial position of a single particle. From the 
definition of the magnetization, once a change of variables has been applied from 
{0,p) to (6*0,^0)7 we obtain: 



and Eq. is recovered. 

The next step is to incorporate the additional contributions, so far neglected. 
In particular, we focus our attention on the terms i^"Ap^"^'^ in Tab. [T] For 
A9 w TT (i.e. Mq <C 1), the dominant term is — 7\/oAp^t^/36 since all the other 
terms are of higher order in Mq. This latter can be seen as originated from a 
modification of Eq. p2)) where a constant factor c is being introduced as: 



The coefficient of ApH^ of Eq. ^ is replaced by ApH^ - 2cAp'^t'^. Therefore 
c = 5/3 matches the dominant term — AfoAp^t*/36, corresponding to the order 
n ~ 2 . The approximation of the magnetization given by Eq. (|14p is in better 
agreement with the numerical simulations. In particular for the position of the 
resonances, Equation ()14|) gives Pr — \/ Ap^ — c, which is closer to numerical 
values as shown Fig. [5] (dashed line). However, this additional term does not 
balance the analogous contributions associated with higher orders (n > 2) for 
which a slightly different value of c is required. Deviations are however reason- 
ably small (less than 10 %) over the range of inspected coefficients. The above 
argument can be extended to the case where AQ < tt, so accounting for the 
terms proportional to Mq: In practice, an additional term of the type ciMq is 
introduced in the square roots of Eq. (|14p where ci is a constant. 

3.2. Low-energy limit 

We now consider the low-energy limit [/ ~ (1 — Mq)/2, that is Ap <C 1. 
We notice that this limit is close to the line which marks the forbidden region 
in the parameter space {Mq,U) (see e.g. [l^). In what follows, we find an 
approximation of the coefficients of M{t) proportional to Ap'^. We first observe 




M{t) = Mosinc(iv/Ap2 _ c). 



(14) 
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that the Liouville operator ([7]) either increases or decreases by one order the 
exponent of p. Thus, the odd powers Ti^M contain a set of elements hn.m with 
m odd. For the waterbag initial distribution, such terms vanish so M (t) is an 
even function. Then, the recursion relation Q is generated by H^. In the 
low-energy limit, if the kinetic terms are neglected, we get 

Ti.'^bufi ~ {bi.obu-1,0 - &-i,o&«+i,o) ■ (15) 

An algebraic expression of the magnetization in the low-energy limit is ob- 
tained by studying the sequence of terms at the lowest order. In this way, we 
approximate M as 

M (t) w J2 (""^1,0 + Pnh-ifih2fl -f 7„6_i,ofo?,o) (16) 

n=0 

Using Eq. (|15p . we deduce that, at a given order the bi.o term comes from 

the bi,Q term at order n, with an+i — Q!„/2. From this recursion relation, we 
deduce the formula 

= ^- (17) 

The 6-1,0^2,0 term at the order n -I- 1 is generated from both the 6i_o and the 
^-1,0^2.0 terms at the lower order n. The recursion relation becomes Pn+i = 
{Pn — Q;„)/2, which leads to 

/3« = -|r- (18) 

The third term in 6-i,o&i o '^'^^ only generated through the reduced operator 
given by Eq. (jlSp . but also from other nonlinearities: The latter terms have been 
neglected in Eq. psp . but appear when considering Ti"* (and possibly higher 
powers of 7i) in the low-energy limit. We resort to an ansatz for ^2m fitting the 
coefficients derived algebraically up to 710: 

for large n. It follows that the magnetization in the low-energy regime is ap- 
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Figure 6: The magnetization M{t) is plotted versus time, in the low-energy regime (U = 0.01), 
for Mq = 0.2 (blue curves) and 0.6 (red curves). The solid lines refer to Eq. I I20I I for kg = 20, 
the dotted ones to direct A''-body simulations. 

proximated by 

M{t) « Mocosh(^-^^ -MoM2(0)^sinh(^-^ 

-Af^icosh y^t^ . (20) 

This expression of the magnetization is compared with numerical simulations in 
Fig. [HI We notice the good agreement up to the saturation regime. As expected, 
for longer times, the approximation gets worse due to higher order nonlinearities 
which have been neglected. 



3.3. Out- J- equilibrium phase transition 

As previously reported, increasing the energy [/ at a fixed value of the initial 
magnetization A/q leads to a drastic change in phase space which materializes 
as an out-of-equilibrium phase transition 4] from an inhomogeneous to a homo- 
geneous phase. This phenomenon was first explained by invoking a principle of 
entropy maximization, based on the so-called theory of "violent relaxation" 



Another dynamical explanation of such transition comes from a bifurcation anal- 
ysis in phase space [sj. 

We here show that the transition can be also retrieved when tracking the 
short-time behavior of the magnetization. It means that the system relaxes 
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very quickly in its metastable phase. The idea goes as follows: We monitor 
the magnetization dynamics via the analytical expression obtained from Eq. ^ 
and store the first local maximum, for each choice of the pair {Mq, U). In case 
the series diverges, without passing through a local maximum, the intensity 
is recorded when its derivative crosses a given threshold (as a polynomial, it 
eventually explodes). We choose fco = 20 in the algebraic computations. The 
resulting values of the magnetizations are displayed in Fig. [7] adopting a color 
code which continuously interpolates between the large (Af w 1) and small 
(M « 0) magnetization. As clearly shown, the upper portion of the parame- 
ter plane corresponds to almost homogeneous configurations while magnetized 
phases are observed as the energy is reduced for fixed AIq. This scenario qual- 
itatively agrees with direct numerical integrations, as confirmed by inspection 
of Fig. [51 In the TV-body simulations, the available parameter space {Mq,U) 
is partitioned in small cells, each associated with a reference water-bag distri- 
bution (that is, a two-level distribution). The QSS magnetization is measured 
by averaging the numerical time series over a finite time window after relax- 
ation. The average QSS magnetization is then represented using the same color 
code as above. When comparing Figs. [7] and [8] it should be emphasized that 
the QSS regime occurs significantly after the violent relaxation process, beyond 
the first local maximum of the magnetization which is computed here. The 
results show that the average magnetization as recorded in the QSS correspond 
approximately to these local maxima in the non-homogeneous phase. A better 
quantitative matching can be obtained by considering higher order terms (larger 
fco). Even though, improving the accuracy of the theoretical analysis is a crucial 
requirement, already at this level of approximation it emerges a phase transition 
as clearly shown in Fig. [71 

A reference line in parameter space (Mq, U) marking the transition between 
the magnetized and unmagnetized phases can be computed based on the cele- 
brated Lynden-Bell procedure, also known as the violent relaxation theory. The 
central idea of the Lynden-Bell approach consists in coarse-graining the micro- 
scopic one-particle distribution function f{0,p; t) by introducing a local average 
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in phase space. Starting from a waterbag initial profile with a uniform distri- 
bution /o , a fermionic-like entropy can be rigorously associated with the coarse 



Xini + (l-X]in|l-i- 



The 



grained profile /, namely s[/] = — JdpdO 
corresponding statistical equilibrium, which applies to the relevant QSS regimes, 
is hence determined by maximizing such an entropy, while imposing the con- 
servation of the Vlasov dynamical invariants: Energy, momentum and norm 
of the distribution. The analysis translates into the out-of-equilibrium phase 
transition line reported in Ref. [4]. Notice that the Lynden-Bell scenario re- 
called above formally applies to the waterbag initial condition from which the 
fermionic principle is derived. Different energy functionals are at variance to 
be assumed when dealing with more complex initial conditions and there is no 
a priori guarantee that the maximum entropy strategy would perform equally 
well. Aiming at extracting a transition line from the viewpoint of the bracket 
calculation, one could impose a critical threshold Mc to the magnetization: First 
local maximum values of the magnetization larger than Mc are assumed to yield 
a magnetized QSS, while for magnetization below the reference value Mc the 
system evolves toward a homogeneous QSS. The (arbitrary) choice Mc = 0.4 
leads to a transition line (dashed line in Fig.[8|) which resembles qualitatively the 
Lynden-Bell line (solid line). Notice that magnetized patches are numerically 
seen to extend over the region of homogeneous QSS, so effectively deforming 
the transition boundary in a non trivial way. Interestingly, such islands are 
entrapped in the wiggles of the bracket transition profile. 

In conclusion, the bracket method returns sensible information on the exis- 
tence of an out-of-equilibrium transition, so resulting in a powerful tool for those 
generalized settings where the Lynden-Bell ansatz proves inadequate (as for in- 
stance, for Gaussian initial conditions) or, at least, cumbersome (e.g. multi-level 
initial distribution). 
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Figure 7: Map of the magnetization evaluated at the first local maximum of Eq. ^ in the 
{Mo,U) plan. The data refer to the theoretical prediction calculated for fco = 20. The white 
region is the forbidden one. 




Figure 8: Map of the QSS magnetization in the (Mq, U) plan, as recorded via direct Af-body 
simulations {N = 10000). The solid line refers to the Lynden-Bell prediction. The dashed line 
stands for the bracket transition line with threshold magnetization set to Mc = 0.4. 
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4. Conclusion 



In this paper we have focused on a paradigmatic Hamiltonian mean-field 
model, often being investigated for its long- living Quasi Stationary States (QSS), 
and for its peculiar out-of-equilibrium dynamics and phase transitions. Starting 
from an out-of-equilibrium initial conditions of the waterbag type, the system 
rapidly evolves toward an intermediate dynamical regime, distinct from the 
corresponding equilibrium configuration. When increasing the number of in- 
teracting elements the time to equilibration gets longer and formally diverges 
when the thermodynamic limit is performed. Aiming at shedding light onto 
the short time dynamics, which ultimately governs the QSS emergence, we have 
here resorted to an analytical approach. The idea is to develop an algebraic 
technique based on the Lie-Poisson structure of the HMF dynamics. In doing 
so we are able to return an analytical prediction for the global magnetization as 
a function of time, a macroscopic parameter sensitive to the microscopic par- 
ticles evolution. Two limiting cases are explicitly considered, respectively the 
high and low energy settings, and shown to yield to tractable expressions for the 
magnetization amount. In general, and due to the perturbative nature of the 
calculation, the full analytic expression contains a vast collection of terms which 
are difficult to handle. The number of terms involved increases rapidly with the 
order of the approximation making it practically difficult to address the dynam- 
ics in the relevant, saturated, QSS regime. However, targeting the analysis to 
the first local maxima in the magnetization, and accounting for 20 orders in the 
perturbative expansion, the existence of an out-of-equilibrium phase transition 
was singled out, separating between homogeneous and non-homogeneous QSS. 
This transition was already recognized in Ref. and interpreted using an ad 
hoc maximum entropy principle suited for waterbag initial profiles. Although 
the calculations are carried out for the so-called waterbag initial condition, the 
technique we use in this article is rather fiexible and can be readily extended 
to other, possibly more general classes of initial conditions so returning fully 
predictive scenarios. We also notice that the proposed method can be adapted 
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to other contexts where long-range many body interactions are at play. The 
method is particularly adapted to short-time dynamics (transients, metastable 
states, violent relaxation, etc.). 
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